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DYNAMICS AND COLUMN DENSITIES OF SMALL PARTICLES 
EJECTED FROM SPACECRAFT 

INTRODUCTION 


Small particles are often generated by a spacecraft in flight and can interfere with 
intended operations in a number of ways. One of the primary areas of concern is the 
light scattered from such particles. Of secondary concern is the molecular material 
sublimating from them. To calculate these effects it is necessary to know the column 
density. 


The column density depends on the source rate and the velocity. To a first order 
approximation the velocity may be taken as the initial velocity of separation between the 
particle and the spacecraft. However, this velocity can be modified significantly by orbital 
dynamic and drag forces. These effects are computed in this study. 


ANALYSIS 


Consider the satellite in a circular orbit with velocity v 0 . Let the particle be 

ejected with a velocity Sv making angle 7 with the velocity vector and T with the 
radius vector. Define the x-axis to be along the velocity vector, the y-axis to be along the 
radius vector, and the z-axis to be opposite the orbital angular momentum. The 
components of 5v are 


->■ 

5v 


|Sv| 


I cos 7 \ 

cos T 

\x / 1 - cos 2 7 - cos 2 T j 


(D 


The particle will have a change in orbital inclination given by 5i = 5v z /v 0 and will 
oscillate about the z direction according to 


Sz = r sin <p 5i 


(2) 


where r is the radius vector and <f> is the angle traveled from the release point. 



By defining a quantity 


e - v 2 /v Q 2 - 1 

where 

v = iv^j + 6v| 



The semimajor axis a and eccentricity e of the particle orbit are 



and 


2 - v 2 /v @ 2 

where 5 is the dispersion angle or 6 = Sv y /v 0 . In terms of e , 
e 2 = 1 -(1 -e 2 )(1 -6 Vy/Ve ) 



(4) 


(5) 


( 6 ) 


2 



and 


Hereafter the subscript e refers to the spacecraft. 

The mean anomaly M 0 at the point of release is obtained by requiring the radius 
vector r to equal r 0 . Using the relation between radius vector and eccentric 

anomaly E and equation (7), 

r @ /a = 1 - e cos E 0 = 1 - e , 


from which 

E 0 = ± cos- 1 (e/e) 


( 8 ) 


where E 0 is the eccentric anomaly at the time of release. 

The sign ambiguity introduced by the arc cos may be resolved by 
taking E 0 positive if the particle is projected in the direction of the radius vector 
(0 < 7 <1 80 deg) and negative otherwise. 

The mean anomaly M 0 is found from the relation 


M 0 = E 0 - e sin E 0 


(9) 


The mean anomaly may be expressed as a function of time by 


M 


M 0 


2 it t 

+ 

T 


( 10 ) 


where T is the period. This period is related to the spacecraft period by 
T - T 0 (a/r = T s (' + f e + • • ) ■ 


( 11 ) 
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The mean anomaly of the spacecraft relative to the particle perigee is just the true 
anomaly v 0 of the particle at the time of release. This is found by requiring 


1 - e 2 r © 

= — = 1 - e 

1 + e cos v 0 a 


( 12 ) 


or 


Vn = ± COS M 


r- ; i 

Led -e)J 


(13) 


The same resolution of the sign ambiguity stated for E 0 applies here. 

The mean anomaly for the spacecraft is given as a function of time as 


~ ‘'o + 


2 tt t 


(14) 


The position of the particle relative to the spacecraft in a space-fixed coordinate system 
with x' in the direction of the particle perigee is given by 


5 x' 
r o 


- e 


[cos M + e (cos 2M - 3) + i e 2 (3 cos 3M - 3 cos M) + . . . j 


- cos M 


e 


(15) 


*£ _ l [ 
^ " 1 - * L 


sin M + i e sin 2M + e 2 (9 sin 3M - 1 5 sin M) + 


] 


- sin M , 
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These may be expressed in a rotating coordinate system such that x is in the direction 
of the spacecraft velocity and y is along the radius vector by rotating the system 
through angle M 0 . 

The drag on the particles is treated by allowing the semimajor axis to vary 
according to 


a 

a 


E 

E 


(16) 


where E is the energy of the orbit. From the Virial Theorem, 


2 < KE ) = - < PE > 


(17) 


where the bracket terms are the time-averaged kinetic and potential energies. Therefore, 
the energy of the orbit is 

E = - t m < v 2 > . < 18) 

The rate of energy loss caused by drag is 


E = - Fv = -C D A p a v 3 ( ,9 ) 

where is the drag coefficient which is the order of unity, A is the frontal area, 
and p a is the atmospheric density. 

For nearly circular orbits the velocity is almost constant so that the difference 
between the average velocity and the rms velocity is negligible. Therefore, 


2 Cp A p a v 


a 


m 


( 20 ) 


It is useful to express the change in a in terms of the traction change per orbit, 
or 


5 



da . 
— = a T 
dn 


where the period T is 


Therefore, 


da 4 * C D A Pa a2 


where n is the revolution number. This has the form 


= - k a 2 


Integrating with the boundary condition a = a 0 at n = 0 , 


I 1 

— - — = k n 
a i a 0 


“ “ 1 + a„ k n • 


Assuming the particles are uniform ice spheres, Cn = 1 , a„ = 6.798 X i 0 6 and 
p a = 4.82 XI O' 1 2 kg/m 3 (at 420 km). 
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4 7r C D A p a a 0 


6.15 X lO' 7 

d (m) 


(27) 


a 0 k 


m 


where d(m) is the particle diameter in meters. Finally, equation (15) becomes 


5 x' 1 - a 0 k n [ 
- cos 

r a 1 “ e 0 L 


1 


cos M + — e (cos 2M - 3) 


+ - e 2 (3 cos 3M - 3 cos M) + 
8 


. . - cos M 0 (28) 


5 y' 1 - a Q k n f 

V 1 - e o L 


sin M + - e sin 2M 

2 


+ — e 2 (9 sin 3M - 1 5 sin M) + 

24 


. . . J - sin M 0 


where a 0 k is given by equation (27), and n is the time in orbital periods. 

To find the distance the particle has moved in one orbit, replace M by M + 2 it 

3 

and M 0 by M 0 + 2 tt T/T 0 = M 0 + 2 tt (1 + - e + ...) in equation (15) and subtract the 
original equation. This yields 


6 x' = a[cos M 0 (1 - cos 3 rr e) - sin M 0 sin 3 it e] 

(29) 

6 x' = a[sin M 0 (i - cos 3 n e) - cos M 0 sin 3 7r e] . 

Since e « 1, (1 -cos37re) involves only second order terms which may be neglected. The 
distance traveled in one orbit AS becomes 


— = (5 x' 2 + 8 y' 2 ) « | a sin 3 7r e I 1 3 it e a | . (30) 

An 
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Under the influence of drag, a changes by a = a 0 (I - a 0 kn) and e changes by 
e = e o ■ a o kn. Both e 0 and a 0 kn « 1 . Ignoring terms of higher than first order, 


AS 

An 


I 3 ir a 0 (e 0 - a 0 k n) | 


(31) 


Integrating, the average separation as a function of time is given by 


S 


3 7 r a 0 


e t a 0 k t 2 
T 2 T 1 


(32) 


where 


a o 


r 


o 


1 - e 0 


eo 


2 


5 v 

— cos 7 


+ 



and T is the orbital period of the particles. Actually, this changes also with a , but the 
difference is only in terms higher than first order. 


RESULTS 


Typical trajectories in the rotating coordinate system are shown in Figures 1 
through 3. Figure 1 shows particle trajectories for an ejection velocity 
of 5v = 0.001 v 0 at various angles and no drag. It is interesting to note that particles 

ejected along the velocity vector pass above the spacecraft and move to the rear, while 
the opposite is true for those ejected to the rear. This can be easily understood from the 
fact that addition of velocity raises the orbit and increases the period; therefore, the 
spacecraft will pass under and ahead of the particle. 


Particles ejected at 90 deg and 270 deg have approximately the same period and 
will, therefore, almost reencounter the spacecraft one orbit later. Also, note the fact that 
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trajectories have a mirror symmetry; i.e., the 225-deg trajectory is the mirror image of 
the 45-deg trajectory. 

Figure 2 shows a family of trajectories for an ejection velocity 
of 5v = 0.0005 v @ and no drag. Trajectories for angles 1 80 deg < 7 <360 deg were 

deleted for clarity since they are mirror images of those shown. The epicycles can be seen 
for the 7 = 60 deg and 1 20 deg case. Figure 3 is the same case as shown in Figure 2 
except that the scale is smaller. The tendency for the particles to spread out along the 
orbit is seen. In the direction of the velocity, the particles tend to stay below the flight 
path, while the opposite is true in the negative velocity direction. 

Figure 4 is a plot of the radial distance between the spacecraft and the particle as 
a function of time. For comparison, the trajectory corresponding to initial ejection 
velocity is shown. It may be seen that particles leave with this velocity and that clearing 
rates may either be enhanced or retarded by the orbital dynamics. The 45- to 225-deg 
trajectory provides the fastest initial separation. A similar set of curves plotted over 
several orbits is shown in Figure 5. 

The effect of drag is the decrease of the semimajor axis. This will ultimately cause 
the particle to move ahead of and below the spacecraft. Figure 6 shows the trajectories 
of particles ranging from 1 00 /am to 1 cm ejected at 7 = 60 deg with a velocity 
of 5v = 0.001 v 0 . The 100-pm particle turns and is rapidly swept away. The 1-mm 

particle first proceeds opposite to the velocity, makes several loops, and then rapidly 
moves in front of and away from the spacecraft. The 1-cm particle moves similarly to the 
drag-free case. It is being retarded and will ultimately reverse direction and pass below 
and in front of the spacecraft. 

Figure 7 is identical to Figure 4 except that drag for a 100-pm particle is 
considered. Now the symmetry is destroyed and the 225-deg ejection angle results in very 
rapid clearing. The 135-deg case has a slower initial clearing rate and the 45- and 315-deg 
cases delay the rapid acceleration from the drag until almost one orbital period. 
Therefore, if overboard dumps are required, there is some advantage to projecting the 
material toward the rear and down. 

Figure 8 shows the separation rates for 100-pm, 1-mm, and 1-cm particles for 
different values of 7 . Note that clearing is always more rapid for particles thrown to the 
rear of the spacecraft. 


COLUMN DENSITIES 


The ultimate use of this study is to compute densities in a column along some 
line of sight. To do this properly, the contribution from a given source must be obtained 
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by following the trajectory of each particle and the probability of the particle being 
found in the specified column obtained. This can be done by Monte Carlo techniques. 
However, for more general estimates simplified models will be used. 

First assume a free expansion model in which particles are emitted uniformly 
from the surface of a sphere with radius R 0 and radial velocity 5v r . The number 
density is just 


N(R) 


Nj 

4 ir R 0 2 S v r 


(33) 


where is the total emission rate (particles per unit time). The column density is 
found by integrating along a column, 


n 


c 


/ N(R)dR = 
Ro 


Nj 


4 ?r R 0 5 v r 


(34) 


This model applies if the 5v r is large enough so that the particles disperse to an 

undetectable concentration in less than one orbit. However, because of the orbital 
motion, the particles do not expand indefinitely but are confined to a tube about the 
orbit. 


For an order of magnitude estimate, assume that all the particles leaving the 
spacecraft are uniformly distributed in cross sector in a circular tube with radius £ and 
whose length is increasing with time at the rate S along the orbit. 

The y-component of particle oscillation is given by 




ae 


r 



(35) 


and the tangential amplitude is 


5 v„ 


( 36 ) 
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The resultant amplitude is 


t = <v + Sz 2 ) 14 


— (4 Sv x 2 * SVy 1 + 6 « z ! )' /! 


(37) 


or 


£ = r., — (3 cos y + 1 ) 1/2 
v ® 


Note that £ varies only by a factor of 2 over the entire range of y . For the 
accuracy required here, it is sufficient to ignore the y dependence and set 


8 v 



(38) 


The number density contribution from particles inserted at angle y anywhere in 
the tube is given by the number inserted during the interval At divided by the volume of 
a length of the tube equal to SAt , or 


A N ( 7 ) _ N ( 7 ) 

7 T £ 2 S At £ 2 S 


where S is given by 


1 AS _ 3 7 T a 0 
T An T 


In terms of S , this becomes 



S = ,/lSo 2 ±2S 0 S | 


(39) 


(40) 


(41) 


11 



where 


Sq 


3 ir a 0 e 0 

T 


and 


So 


3 7 r a 0 2 k 
T 2 


The positive sign which decreases the column density applies to particles moving ahead of 
the spacecraft, i.e., those ejected to the rear with 7 > ir/2 resulting in e < 0 . 

The column density may be calculated by integrating ? 7 (S) along the sight vector. 
Let the sight vector make angle \p with the orbital path. The contribution from particles 
released with angle y to the column density is 


S, 


r ds 

5 n c = / „(S) 

c Q cos \p 


N(T) 


S, 


f ds 

■I 7s 2 ±?.s 


jr ? cos $ J Q [S 0 2 ± 2 S 0 S] l /2 


(42) 


S 0 may be expressed in terms of y . From equations (30) and (3), 


• 37 ra 0 / 5v\5v 

S 0 = -r- (2 cos 7 + — 


v o/ V 0 


(43) 


or, using the fact that T = 27rr 0 /v @ , 


3 / 5 \ 

S 0 = - - S v 2 cos 7 + — 

2 \ v o/ 


(44) 


The sign was lost in the derivation when extracting the root; however, it is clear 
that particles ejected with 7 <ir /2 will move behind the spacecraft, whereas 7 > 7 r /2 
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results in positive velocities. Some care is required to insure the proper sign in equation 
(42). For \p<nl2 , particles ejected with y>n/2 will be seen directly, and they will be 
accelerated by drag. Therefore for < n/2 , 


n, 


rr/2 ^ S, 

^ N ( 7 ) cry ^ ds 

0 7T p cos 4, 5 [S 0 2 + 2 S 0 si 1/2 


(45) 


where S t is the distance at which the line of sight intersects the tube containing the 
particles. This is 


S, 


£ cot \}/ 


whichever less. 


(46) 




Particles ejected with y<n/2 will initially move opposite to the velocity vector. 
They will be retarded by drag and will ultimately turn and be accelerated along the 
velocity vector. When they pass the spacecraft, from equation (41), it may be seen that 
theii; separation velocity will be equal in magnitude to S 0 , but they will be moving in the 
opposite direction. From this point on they will behave just like particles ejected 
at jt/2 - 7 . Therefore, the total column density for \p < tt/ 2 will be twice the integral of 
equation (45) over 0 < 7 < ir/2. 

For i//> 7t/2 , the situation is somewhat more complicated. There will be no 
contribution from particles ejected with 7 > tt/2 since they will drift in front of the 
spacecraft . 1 For y<w/2 , the particles will drift to the rear. The effect of drag is initially 
to slow the particles, which requires that the terms in the denominator are subtractive. 
This will come about naturally since S will be negative. The limit Sj will be that 
specified in equation (46) with opposite sign, or the turning point of the particle, 
Sj = -S 0 2 /2S 0 , whichever occurs first. The returning particles will give a contribution 
equal to those going away. Therefore, the total will be twice the amount found by 
integrating equation (45) from 0 < 7 < tt/2. 


1. Particles ejected with y>n/2 will, of course, initially pass through the field of view 
as they make their initial loop around the spacecraft, but the long-term motion is being 

considered here. 
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Integrating equation (45) once yields 


n 


c 


$ 2 N (7) d7 _2 
i it? cos 1 p S 0 


[(So 2 + So S ,) 1/2 



(47) 


l : or : p<n/2, the above holds for all 7 . If the spacecraft is emitting particles isotropically 
at a rate N t /4tt particles/sec/steradian, the number density can be found by 
replacing N( 7 ) by N T dft/4jr or ViUj sin7 , multiplying by 2 and integrating 

from 0 < 7 < tt/ 2 . This results in 


3 Nj 6v 
4 ir £ 2 cos 4* S 0 


(2 + 5v/v 0 ) 2 (A, - 1) + 



CA 2 - 1) 


4 S 0 S] 
+ 9 5v 2 


In 


"(2 + 5v/v 0 )(l + A,)' 


5v/v 0 (1 + A 2 ) 

J 


(48) 


where 


4 S 0 S, 

9 5v 2 (2 + 6v/v 0 ) 2 


4 S 0 S, 

1 + 9 5v 2 Sv 2 /v 0 2 


For ip > 7r/2 , the above holds for particles that do not turn before - Sj , or for 

which 



£ cos 1 p 


whichever greater. 


(49) 
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For particles that turn before they move out of the viewing column, Sj must be replaced 
by -S 0 2 /2S 0 in equation (47). This applies to all particles for which 


S 0 2 < 2 S 0 £ I cot yp | 


(50) 


or for which 


r 


o 


k > 


3 it [2 cos 7 5 v/v 0 + (5v/v 0 ) 2 ] 2 
2 Sv/v 0 | cot ^ | 


(51) 


At Skylab heights (420 km) the drag term r 0 k for a spherical particle with 
density of 1 gm/cm 3 is 0.615/d (nm) where dfcmi) is the diameter in microns. This 
means that all particles smaller than 


d(/um) < 


0.615 2 5v/v 0 1 cot yp | 

3 7 r [2 Sv/v 0 cos 7 + (Sv/v 0 ) 2 ] 


(52) 


will turn before they move out of the viewing column. For 7 -O and 
Sv/v 0 = 0.001 (7.8 m/sec). 


d(pm) < 32.6 j cot ^ | 


Therefore, integrating equation (47) with S, =-S 0 2 /2S 0 applies to all cases for which 
the viewing angle is greater than \p > cot 1 [-d0rm)/32.6] . For 100-fzm particles, this 
restricts to greater than 162 deg. For 10-pm particles, ^ m i n is only restricted to 

greater than 107 deg. Integrating equation (47) for yp > ^ m j n and doubling to account 

for the particles after they turn yields 

3 y/2 N t 5v 
c 2 7 T £ 2 cos yp S 0 


5v 

. + — 
v, 


1 ] 


(53) 
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The intermediate cases, 7 t/ 2 <! must be treated by integrating equation 

(47) from 0 to 7 , with S, given by £coti p and from 7 , to ir/2 with S, = -S 0 2 /2S 0 
where 7 t is the value for which equation (52) is an equality. 


The column densities for a spacecraft emitting particles isotropically are plotted as 
a function of velocity in Figures 9, 10, and 11 for look angles \p = 0 , 7 r /2 , and n. The 
column densities are normalized by the column density that would result in the absence 
of orbital and drag effects, equation (30). The drag-free case is obtained by taking 
the lim as S 0 -> 0 of equation (47). This yields 


n 


c 


n t s, 

6 jt £ 2 cos \jj 8\ 


In 


2 + Sv/vA 

, 5v / v ® / 


(54) 


For the \jj 0 case, it may be seen that orbital effects on the column density are 
not important unless 5v < 4 m/sec. For small particles which produce the highest 
scattering efficiency, drag effects rapidly clear the particles even if their initial velocity is 
low. The fact that at high ejection velocities the column density from small particles 
exceeds the no-drag case is because of the factor of two resulting from the drag 
eventually causing all particles to move into the field of view. 

For the \}/ = ir/2 case, the orbital effects are negligible for all particles with 
velocities more than a few tens of centimeters per second. With drag effects, the smaller 
particles are accelerated out of the field of view so rapidly that they do not contribute to 
the column density regardless of their initial velocity. 

For the direction 4/ = ir the effect of drag is to greatly increase the column 
density, particularly for the larger particles. This is because of the increased number 
density brought about by the particle’s low velocity in the vicinity of the turning point 
Equation (53) is restricted to values of 5v small enough to cause all particles to turn 

before s/2^ . These velocities are indicated by the termination of the plotted lines. 
As Sv is made large, these curves will approach a value of twice the drag-free case. 


Finally, Figure 12 shows the relative column density for the Apollo in lunar orbit. 
It was assumed that the Apollo was an isotropic source with an effective radius of 2 m 
Shown for comparison is the tfr = 0 , 90 deg no-drag cases for Skylab in a 420-km earth 

orbit. Note that in lunar orbit, orbital effects are not significant for ejection velocities 
more than 0.56 m/sec. 
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SUMMARY 


The effects of orbital dynamics and drag on small particles released from 
spacecraft in circular orbits have been calculated. It is shown that such particles will 
become distributed along the orbital path and that the net separation rate between the 
particle and the spacecraft is given by 


S 



5v 

2 — cos 7 + 



2 C D A p a v 0 t 
m 


where v 0 is orbital velocity, 8v is ejection velocity, y is the angle of ejection relative 
to the orbital velocity vector, Cp is the drag coefficient, A/m is the area to mass ratio 
of the particle, p a is the atmospheric density, and t is the time after ejection. Particles 

ejected randomly will be confined by orbital dynamics to a tube of cross sectional area 
given approximately by ir£ 2 where £ = r 0 5v/v 0 and r 0 is the orbit radius vector. 

The column densities resulting from a continuous isotropic release of particles in 
circular orbit from the Skylab were compared to a similar release in the absence of 
orbital and drag effects. It was found that substantial increases in column density along 
the antivelocity vector can result from orbital dynamics and drag forces for particles 
larger than 10 pm ejected at velocities less than 4 m/sec. Drag and orbital effects 
significantly reduce the column density when viewing at 90 deg to the velocity vector. 
Drag forces produce rapid clearing of particles smaller than 100 pm in front of the 
spacecraft and subsequently reduce the column density when viewing along the velocity 
vector. Larger particles ejected at velocities less than 1 m/sec can cause some increase in 
the column density because of their longer time in the field of view resulting from their 
confinement along the orbital path. 

Finally, the case for an Apollo spacecraft orbiting the moon was considered. 
Despite the lack of drag, there is negligible accumulauon of particles even in a column 
along the velocity vector for ejection velocities larger than 0.5 m/sec. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama, September 1973 
502-21-28-0000 
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les at 5v = 0.001 v 0 for no drag. 
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Figure 3. Particle trajectories of Figure 2 with reduced scale to show the relative motion 

of the particles along the orbit. 
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Figure 4. Separation distance between particle and spacecraft. (Initially the particles leave at their given separation 
velocity, but after ~ 0.1 orbit the orbital dynamics may either accelerate or retard the 
separation rate, depending on the ejection angle.) 
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Figure 5. Separation distance between particle and spacecraft shown in Figure 4 but with 

reduced scale to show long-term behavior. 
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Figure 6 . Particle trajectories for 7 = 60 deg, Sv = 0.001 v 0 , with drag at 420 km. (The 100-jum particle turns and is 

rapidly swept away. The 1-mm particle moves initially to the rear of the spacecraft, turns after several orbits, and 
then moves in front and separates rapidly. The 1 -cm trajectory is similar to the drag-free case, but the effect 
of drag is seen to retard the separation. This particle will also eventually turn and move 

in front of the spacecraft.) 
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Figure 8. Separation distances for 1-^m, 1-mm, and 1-cm part 
that clearing time is decreased by projecti 
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Figure 8. (Concluded) 
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Figure 9. Column densities from particles released continuously and isotropically with velocity 5v in a 420-km circular 
orbit looking along the velocity vector compared with column densities that would result from a free expansion 
of the same particles. (There is no significant accumulation from orbital effects for release velocities larger 
than 4 m/sec. Smaller particles are rapidly cleared by drag and particle sizes below 100 jum will not 

accumulate significantly for velocities above 18 cm/sec.) 
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Figure 1 1. Similar to Figure 9 except that the viewing direction is opposite to the velocity vector. (In this case substantial 
accumulation can take place because the particles moving rearward are retarded by drag and eventually accelerated 
back toward the spacecraft. The lines are terminated where equation (53) no longer applies because particles 
with larger Sv will move out of the field of view before they turn. The column densities for larger 5v 

will eventually approach twice the drag-free value.) 
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Figure 12. Column densities from orbital effects (no drag) in earth orbit (S 

densities from the same particle release in lunar orbit 





